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Abstract 

We demonstrate that the main part of the observed nuclear suppression of dilepton 
and charmonium production, especially the x-dependence of the suppression, can be 
explained by the absorption of fast initial partons. We assume the absorption to be 
xi-dependent and determine the empirical form of this dependence. Several reactions 
at different energies are described with the same parametrization of the absorption 
cross section. The factorization theorem constraints and some alternative models are 
discussed. The obtained description of data allows to make conclusions about the 
origin of nuclear effects in hadro- and electroproduction. In particular, it is concluded 
that the parton recombination and the excess pion contribution are unimportant in 
nuclei and an indication of excess gluons at X2 ~ 0.1 is observed. 

PACS numbers: 24.85. +p, 25.75.+r, 12.38.Aw, 13.85.Ni 
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I. INTRODUCTION 

Several effects have been considered in order to explain the observed strong suppression 
of J/if>— production in hadron-nucleus collisions [|. We can mention some of them: the 
intrinsic charm[Q, the comover interaction 0], the recombination of nuclear partons |], [J, the 
interaction of non-singlet final parton configurations 0, the elastic rescattering of initial 
partons 0, the energy loss of initial partons |J and the absorption of initial partons | ID]. 



However, we agree with the authors of Ref.JTTJ that the x— dependence of charmonium 



production suppression has not yet been described. Without a reliable description of the 
charm production suppression in hadron-nucleus interactions, it is hard to analyze this effect 
in heavy ion collisions while searching for the quark-gluon plasma signals [|12[]. Moreover, 
a nuclear suppression has been observed in the Drell-Yan (DY) production by pionsfTB] 



and, though much less, by protons [fL4j. A self-consistent description of hadroproduction 



suppression in nuclei, as well as of the nuclear effects in electroproduction|T5[, is now one 
of the important goals of high energy nuclear physics. We here develop a model, which 
successfully describes several reactions at different energies and, as we try to argue, does not 
contradict the true factorization forecasts. 

We now briefly overview the status of the above mentioned models. 1) The intrinsic 
charm model is, in principle, able to account for the observed xp— dependence of the J/ip— 
suppression, namely the suppression increasing with xp. However, the experimental searches 



for the intrinsic charm in protons[16[ have not revealed any charm signals. With the exper- 
imental upper limit on the amount of the intrinsic charm in protons this mechanism 
cannot explain the observed charm production suppression. Moreover, this mechanism is 
irrelevant for the DY production by pion, where the suppression is also significant, as it will 
be shown below. 2) The comover interaction is most effective at negative xp , because the 
density of comovers increases for less energetic produced quarkonia. In fact, a strong sup- 
pression was observed in quarkonium production at negative xjr ||17|| . In the present article we 
are concerned with the region of positive medium and large xp of produced charmonia, where 
the comover interaction must be minor. 3) In agreement with the strong factorization, the 
parton recombination models suggest that the initial state effects in hadroproduction are in- 



eluded in nuclear structure functions. It has been assumed that nuclear structure functions 
are suppressed at small x 2 due to QCD recombination of partons from neighbour nucleons. 
This mechanism has been widely used to explain the nuclear shadowing in electroproduction 
at low x 2 . This effect would lead to an x 2 — scaled nuclear suppression of quarkonium and 
dilepton hadroproduction. However, a comparative analysis of the data at different proton 
energies clearly shows [TT| that the J/4>— suppression is scaled on xf or x%, but not on 



x 2 . It seems that these data create serious troubles for the nuclear parton recombination 
hypothesis, as we will argue below. 4) Final state interactions of produced non-singlet par- 
ton configurations can be responsible for a suppression of quarkonium production, because 
these interactions are not restricted by the color transparency. However, there is no physical 
reason why this effect should increase with xp at large xp. Moreover, this effect cannot 
explain the DY production suppression. Therefore the final state interactions of produced 
heavy quarks cannot be the only mechanism of nuclear suppression. 5) The elastic rescatter- 
ing of initial partons naturally explains the dilepton and quarkonium yield increasing with 
pp. However this mechanism alone, if not completed by the energy loss of partons, makes 
no predictions for the x— dependence of hadroproduction suppression. 6) The two remaining 
models, assuming the energy loss (EL) of fast initial partons due to soft elastic rescatteringU 
and the absorption of those partons|H|, will be considered below in more detail. We first 



discuss these two models with respect to the factorization theorem. 

The cross section factorization in large-Q 2 reactions can be represented, in the case of 



DY production, in the following form[18, O 



da DY ~ /' ^^f? Y (xi) H(x u x 2 ,Q>) f 2 DY (x 2 ), (1) 

JO X\ x 2 



where f DY (x) are parton densities containing all non-perturbative long-distance effects and 
H, being a projectile- and target-independent function, describes the hard parton-parton 
subprocess. The integration over transverse momenta and the Q 2 — evolution of structure 
functions are implied in Eq.(l). In the lab. frame, the subscripts 1 or 2 correspond to 
projectile or target, respectively. In a general case, when the functions f(x) depend on the 
reaction under consideration, the form (1) corresponds to the weak factorization. The strong 



factorization means that the structure functions are reaction- independent, e.g. the structure 
functions measured in electroproduction are equal to f DY {x) (integrated over transverse 
momenta). The limitations of the factorization have been discussed in several papers (see, 
e.g., Refs.[|19|1-|[22|). In particular, it was obtained^, that in the case of collinear active- 
spectator parton interactions, followed by the productive annihilation of the active parton, 
the strong factorization is valid if 

Q 2 > x 2 LMl 2 T , (2) 

where L and M are the target length (in the rest frame) and mass and It is the transverse 
momentum, transferred in the active-spectator interaction. From Eq.(2), it is reasonable to 
expect that the corrections to factorization, in the case of collinear active-spectator rescat- 



tering, should vanish as 1/Q 2 . It was also concluded |22| that the maximum energy loss of 
an incident parton due to induced gluon radiation, which is expected to be the main source 
of energy losses, is 

A maa; xi = k T L/2Ei ab , (3) 
where is the transverse momentum of the radiated gluon and E^ is the projectile energy. 



This corresponds to an average parton energy loss per unit length dE/dz ~ 0.3GeV/ ,fmp2 
The limitations (2) and (3) directly concern the EL-model, because this model assumes the 
initial state interactions of the type studied in Refs. pl|, [19], |22f . This means that in the 
parametrization for the average parton energy loss in nuclei [[J], 

Axt = Kx.QiQo/QrA 1 / 3 (4) 

{C g = 3 for gluons, C q = 4/3 for partons and Qo ~ 4GeV), it is more reasonable to use 
n=2. This choice is also supported by the observation that the leading contribution to the 
DY production is twist 4. The limitation (3) was significantly violated with n ^0.003, used 
in Ref. |J to fit the data for the charmonium suppression. Nevertheless we will present below 
some numerical results of the EL-model, considering this model as a possible alternative for 
our model. 

The absorption of fast initial partons by nuclear nucleons corresponds to all parton-parton 
processes, which remove the active parton from the production channel. These processes are 



not followed by the productive annihilation of the active parton, in contrast to initial state 
interactions assumed in the EL-model. An example of such processes is a particle produc- 
tion on front nucleons: though the formation of produced particles takes a long time, the 
momentum transfer in productive parton-parton collisions can be quite large, of the order of 
the projectile momentum P. The last statement follows from the observation that at high 
energies many of produced particles are concentrated around the central rapidity New 
particles with central rapidity cannot be produced as a result of only soft parton-parton colli- 
sions. After a large momentum loss, an active parton can hardly participate in the production 
on backward target nucleons. In the terminology of Ref.|l9[], the absorption of active partons 



can be described by hard active-spectator interactions or collinear active-spectator interac- 
tions with large momentum transfer l + ~ P. We emphasize that the proofs of the strong 
factorization performed in Refs.[|l8|]-[^I| do not concern these two types of active-spectator 



interactions and the target-length condition (2) was not proved for these initial state inter- 
actions. For example, the strong factorization has not been proved when l + ~ P/ML (see 
Eq.(2.24) from Ref. ||19||). The cancellation of two-step diagrams, demonstrated in Ref . |20f . 
is not effective for these interactions because the coherency of parton-parton interactions in 
nuclei is lost when / + is much larger than the inverse internucleon distance pi], |IOf . Fur- 
thermore, the limitation (3), obtained for the induced gluon radiation, was not proved for 
some other initial state interactions, for example for the elastic active-spectator scattering 
with large l + . However, it is reasonable to expect that the contribution of large momentum 
transfer collisions is small compared to total hadron-hadron cross sections. This limitation 
will be fulfilled in our model: the parton-nucleon absorption we introduce manifests itself in 
an average proton-proton cross section a^ s < Imb. The large effect of the parton absorp- 
tion in nuclear hadroproduction reactions is explained by the assumed Xi-dependence of the 
absorption: it is concentrated at large x\. In some cases the hadroproduction is sensitive to 
large x±, but the contribution of the large-xi region to the total hadron-hadron cross section 
is small because structure functions vanish at large x. Note that the parton absorption we 
discuss has a negligible effect for one-nucleon targets and the strong factorization is fulfilled 
in this case, in agreement with well known experimental results. On the other side, the 



violation of strong factorization in some hadron- nucleus reactions is an experimental fact, as 
it will be shown below. We will introduce a physically motivated parametrization to fit this 
violation in one reaction with a better statistics and more sensitive to the initial state ab- 
sorption (J/ip— production by 800 GeV protons) in order to study nuclear effects in all other 
reactions. In other words, we will assume that not the strong but the weak factorization is 
valid for nuclear targets. 
II. THE MODEL 

As we argued above, the initial state parton-nucleon absorption can take place due to 
hard and semi-hard parton-parton interactions. Therefore it is reasonable to assume that 
the inclusive parton-nucleon absorption cross section is a function of x\, i.e. a l abs = a l abs (xi), 
where % = g,q denotes gluons or quarks. The weak factorization of hadron-nucleus cross 
sections can now be reformulated as (see Eq.(l)) 

f 1 DY (xi)=F DY (x 1 ,A)f 1 (x 1 ), (5) 

where fi{xi) is the "intrinsic" projectile parton distribution and F is the factor describ- 
ing the parton absorption effects. In quarkonium production, the factor F® can be also 
xp— dependent due to final state interactions of produced partons. The final state interac- 
tions can account for possible comover interactions Q and/or interactions of non-singlet con- 
figurations 0. However, as we mentioned above, these final state interactions are expected 
to change slowly in the positive-Xi? region we are concerned with. Therefore we will assume 
that an x— dependence of suppression factors F is accounted for by the x\— dependence of 
a l abs (xi). The final state effects in quarkonium production will be represented by a constant 
cross section o$. In our earlier workflTOfl, we have used constant cross sections <y l ahs to de- 
scribe initial state effects and the results of the present paper will significantly differ from 



the results of Ref. 10 



In the case of the square well form of nuclear density, the factor F® for the quarkonium 
production (<jf ^ 0) is given by the eikonal form 

F?(x 1 ) = ^--{^^[1 - e-&{l + c£j] - J[l - e-°/(l + C/ )]}, (6) 

°7 ~~ a abs\ X i) C abs Cf 



where c a b s j = 2pa a b s jRA and the nuclear density and radius are denoted by p and Ra- For 
the DY production (a/ = 0), this equation can be rewritten as 

F? ¥ (xi) = i[^f + e- c -(l + O " !]■ (7) 

In Eqs.(6) and (7) we assumed that the exclusive production cross section is much smaller 
than al bs or a f . 

We have to determine the functions a l abs (xi). The physical hadron-hadron absorption 
cross section is given by 

E Jdx 1 a: bs (x 1 )nx 1 ). (8) 

i=g,u,d... 

Since parton distributions diverge as rr — > 0, we must assume a l abs (x\ = 0) = in order to 
obtain a finite cr^. A phenomenological analysis of the J ftp— suppression has shown that 
the data are better discribed when da l abs (xi = 0)/dx 1 = as well. In this case we must 
assume that da l abs (xi = l)/dx\ = 0. This condition can be obtained from Eq.(8) using 
the probabilistic interpretation of parton distributions f l (x). Taking into account these 
limitations, we found that the functions 

< bs ( x i) = C™^^). (9) 

fit well the empirical xf— dependence of the J/tp— suppression in proton-nucleus collisions 
at 800 GeV. We will use the form (9) to derscribe all reactions. The function sin 6 (x^) 
is shown in Fig.l. The parameters o l max characterize the strength of absorption and are 
extreme values of hadron-hadron absorption cross sections when the projectile hadron is 
composed of only one parton i (when x± — 1 for this parton). It is clear from Eq.(8) and 
Fig.l that the physical value a^ bs is much smaller than a l max {^ s ja % max ~ 2%). We may 
assume cr-^. = \o 9 max for all light quark flavors because of the standard color factors. We 
now have to fix the projectile dependence of o % max . At high energies, a pion-hadron cross 
section is about one half of the corresponding proton-hadron cross section. Since <7 % max have 
the meaning of extreme values of hadron-hadron cross sections, we extend the above relation 
between pionic and nucleonic cross sections on a l max as well. For definiteness, the projectile 
dependence will be represented by a l max {pions) = 0Aa l max (protons). As a result, we have 
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only one free parameter, say a^ nax (protons), to describe the initial state effects. We will use 

a 9 max {protons) = 30mb. (10) 

The remaining parameters (a^^protons) and cr%f ax (pions)) can now be calculated. The 
parametrization (9) and (10) will not be adjusted for different projectile energies or for 
different reactions. 

For the quarkonium production, we have to take into account the final state interactions. 
In the study of the empirical regularities in the nuclear charmonium suppression ||11[| , it was 
found that the x— independent part of nuclear suppression can be accounted for by <7/ ~ 
5.8 mb. In our model, these value will simulate the x— independent final state interaction of 
produced heavy quarks. A self-consistent description of charmonium and dilepton production 
will prove a posteriori that in our model there is no room for a significant xf— dependence 
of o$ in the region xp > 0.2. Final state effects have indeed a very modest xp— dependence 
in some models (see, e.g., Ref.|J). The fact that a quite large value of <jf fits the final state 
interactions may indicate that heavy quark-antiquark pairs pass through the nucleus as non- 
singlet configurations, in agreement with the conjecture of Ref.|F|. In fact, o~f is almost three 
times larger than the geometrical-size estimation for the J/ip— proton total cross section p4 



A similar result for the final state cc— nucleon cross section has been obtained in some other 
models||. We think that the symmetrical description of initial and final state interactions 
in quarkonium production, adopted in Ref. 0, is not justified. There can be an important 
difference between the projectile hadron and the produced heavy quark configuration, which 
is not yet hadronized during its propagation in nuclear matter. As we will show, the price 
for the charm suppression fit of Ref. is the underestimation by that model of initial state 
effects in the DY production by pions. 

III. CHARMONIUM PRODUCTION 

To calculate the charmonium production cross section, we take into account both the 
gluon-gluon fusion and quark-antiquark annihilation subprocesses. From the numerical re- 
sults, we concluded that both channels are important in the considered xf— region and the 
relative contribution of these two channels is in agreement with the result of Ref.||. The 
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proton-nucleus cross section is given by 



dQHx F 9 Q 2 {x l +x 2 y Ly 1/yzv z/ 9 Q 2 (^i+a;2) / tr d 

(11) 



The hard partonic cross sections <r M are taken from Ref.[25|. It is known that the nu- 
clear structure functions, ^2(^2); 92(^2) an d <?2( x 2), are not equal to free nucleon structure 
functions. Nuclear structure functions measured in deep inelastic scattering have two pro- 
nounced pecularities[|l5|]: the shadowing at very small x 2 and the depletion (the EMC-effect) 
at £2 >0.4 (the antishadowing at x 2 ~ 0.1 is very small and will not be taken into account). 
The similar effects may be expected in the nuclear gluon distributions. In this section, we 
will not take into account these nuclear effects for the following reasons. If the nuclear shad- 
owing at small x 2 is due to projectile initial state interactions in nuclear matter (in the case of 
electroproduction at small x 2 the projectile can be a virtual quark or meson|26|]), then these 
effects are already included in factors F. We will argue that the alternative interpretation 
of nuclear shadowing, assuming the parton recombination in nuclei, is inconsistent with the 
charmonium production data. The EMC-effect is important at large x 2 , far from the region 
x 2 < 0.14 probed by the current experiments]!], |2[. The latter argument is not correct for 
the quarkonium production in the region x 2 ~ 0.3 studied in some recent experiments [p~7 1 . 
At such x 2 the EMC-effect in nuclear gluon and quark distributions can be noticeable. In 
this section we use the free proton structure functions from Ref.|27| for both f 2 (x) and fi(x). 

The x 2 — and xf~ dependence of nuclear effectiveness a (a = ln{aA/o- p ) /ln(A)) for 800 
and 200 GeV protons is shown in Figs. 2 and 3. The agreement with the data at 800 GeV 
(the solid lines and the diamonds) is not surprising because we used this set of data to find 
the parametrization (9) and (10). The same for 200 GeV protons is shown by the short- 
dashed lines and crosses. This result is already non-trivial, espessially the x 2 — dependence 
of a at 200 GeV. The net effect of final state interactions is represented by the solid and 
short-dashed lines in the small xf— region in Fig.3. 

The parton-recombination model prediction is shown by the dotted line in Fig.2. In 
this case, we have calculated a using the quark and gluon nuclear distributions from Ref. M 
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(the set 1 for gluons) and the constant value <jf = 5.8mb to approximately account for 
the final state effects. Note that the dotted line is expected to describe the data for both 
energies, within a kinematically allowed region of x 2 for each energy, because the difference 
between the results for two energies is less than 1% at the same x 2 . We can see from 
Fig. 2 that the parton-recombination model is inconsistent with the data at both energies. 
Moreover, the parton-recombination model is ruled out by the data, if the origin of the 
nuclear suppression and its gross xp— dependence are the same at both proton energies 
(this is very likely given the results presented in Fig. 3). In fact, in the whole X2— region 
measured at 200 GeV, 0.037 < X2 < 0.14, the effect of the parton recombination is negligible. 
The same should be valid at 800 GeV in the region 0.01 < x 2 < 0.04 since these two 
x 2 — regions correspond to the same xp— region, taking into account the proton energy. . 
It follows from Fig. 2 that the dotted line is in conflict with this expectation. Therefore 
the parton recombination model not only underestimates the charmonium suppression but 
also demonstrates a wrong x 2 — dependence of the suppression at different projectile energies. 
This can be corrected only if the parton recombination model will be combined with some 
other model of nuclear suppression, which predicts a specific dependence on the projectile 
energy at fixed xp. For example, the EL-model does not have this property and cannot 
restore the correct x 2 — dependence, if being combined with the parton recombination model. 
We may assume that the nuclear parton recombination is less important than it is usually 
believed and will not take into account this effect in what follows. 

Though our solid and short-dashed lines in Figs. 2 and 3 correctly reproduce the trend 
of the data, our results underestimate a at x 2 > 0.06. This discrepancy is eliminated if we 
assume that there is a 20%-enhancement of gluons in tungsten at 0.06 < x 2 < 0.15. The 
calculated a with excess gluons included is shown by the long-dashed lines in Figs. 2 and 
3. The resulting \ 2 of our calculations is about 1 at both proton energies. Note that this 
excess gluon contribution affects a only at 200 GeV. The gluon enhancement we assumed 
is in qualitative agreement with the electroproduction data[^|. The momentum fraction 
carried by the excess gluons is about 3% for tungsten and is in qualitative agreement with 



the earlier estimations |59]. It may be assumed that the excess gluons play the role the excess 
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pions have been expected to play - they are responsible for the nuclear binding[29]. Observe 
that the distances involved by the excess gluons are of the order of internucleon distances 
in nuclei. It is also remarkable that the momentum fraction carried by excess gluons is 
comparable with the parton momentum fraction, lost because of the nuclear binding j3(J. Is 
it an indication on the gluonic origin of nuclear forces? We will come back to this problem 
in Sec.V. Note that the excess gluons are not the only possible explanation of the above 
discrepancy. For example, a similar effect can be reproduced by an energy dependence of 
final state interactions. 

The projectile dependence of the charmonium suppression is illustrated in Fig.4. The 
xf~ dependence of a for 200 GeV 7r~— mesons without excess gluons (long-dashed line) is 
compared to the data for pions at this energy (diamonds). The same for protons at 200 GeV 
is shown by the short-dashed line and by crosses. It follows from this figure that our model 
correctly reproduces the projectile dependence of nuclear suppression. Remember that in 
our model, the projectile dependence is represented by the ratio a l max (pions) / a l abs (protons) , 
introduced in Sec. II. However, the long-dashed line still overestimates the suppression. The 
agreement with the data is improved when the excess gluon contribution is included (the 
solid line) in the same way as for proton-nucleus collisions. 

To summarize, the absorption of initial partons can provide a self-consistent description of 
X2~ and xp— dependences of nuclear charm suppression at different energies and for different 
projectiles. In the region 0.06 < X2 < 0.15, the agreement with the data is improved by the 
inclusion of excess gluons. 

IV. DY PRODUCTION BY PIONS 

In the dilepton production by pions, the dominant contribution is due to quark-antiquark 
annihilation subprocess. Since projectile pions contain valence quarks and antiquarks, this 
reaction probes essentially the same target structure functions as the deep inelastic lepton 
scattering. In the case of exact strong factorization, the ratios R of nuclear to nucleon cross 
sections should coincide for these two reactions. Therefore the DY production by pions 
is an excelent test of strong factorization for nuclear targets. We now show that the weak 
factorization is more consistent with the data than the strong factorization. The cross section 
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for producing dilepton pairs in collisions of pions and nuclei is given by 

The suppression factor F DY , given by Eq.(7), was calculated as described in Sec. II. The 



data for this reaction [JTjJ covers a large region of Xi and X2 for two pion energies, 140 and 
286 GeV. Therefore we now have to take into account the nuclear modification of structure 
functions q£ (#2) and q{ (^2)- This modification will be represented by the depletion of nuclear 
structure functions at medium x 2 (the EMC-effect). The form of this depletion that we 



used in numerical calculations reproduces the electroproduction data for nuclear targets |0| . 
As it was discussed above, in the present model there is no physical reason to take into 
account the shadowing of nuclear structure functions at very small X2 since such effects are 
assumed to be included in F DY . We used the free nucleon structure functions from Ref . ||27 1 . 
In numerical calculations, we integrated over the Q 2 — and xp— regions, measured in the 



current experiment p~3 



The X2— dependence of the ratio R for the tungsten at 140 and 286 GeV is shown in Figs. 
5 and 6, respectively. Our results are represented by the solid lines, the predictions of the 
EL-model by the long-dashed lines and the strong factorization result (the EMC-ratio) by 
the short-dashed lines. Here and hereafter the EL-model predictions are calculated using 
Eq.(4) with k = 0.003 and n — 2. The strong factorization is completely inconsistent with 
the data at 140 GeV {\ 2 = 5). The EL-model (\ 2 = 2.3 at 140 GeV) could explain the data 
with a larger parameter k, but this would lead to a futher violation of the limitation (3). 
Our results are in a reasonable accord with the data at 140 GeV (x 2 = 1.1). 

The Fig. 6 is less conclusive because the error bars in the most important region X2 > 0.3 
are larger than the difference between the theoretical predictions and \ 2 < 2 for all models. 
By comparing Figs. 5 and 6, we conclude that in this reaction it is confusing to combine 
the data points for different energies because the X2— dependence of R for two energies is 
very different. For example in the region 0.1 < X2 < 0.3, the suppression at 140 GeV is 
much stronger than at 286 GeV. The "diffraction" minimum at X2 ~ 0.37 in our predictions 
for E n = 286 GeV takes place because the lower mesured mass region, (A.2GeV) 2 < Q 2 < 
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(8.5GeV) 2 , is becoming kinematically forbidden at that Xi- Observe that our model predicts 
a fast depletion at small £2 because larger values of x\ contribute in this region. New data 
at £2 < 0.1 could help to test our model. As it will be shown in the next section, such 
depletion at small £2 has been observed in the DY production by protons. 

The X\— dependence of R is shown in Figs. 7 and 8. Taking into account the very large 
uncertainty in determining x\ for the data points (not represented in these figures but shown 
in Ref.[l3|]), it is hard to make a preference between the models. A better statistics is needed 
to make a definite conclusion about the X\— dependence of R. 

The dependence of R on the mass of the produced dilepton pair is shown in Figs. 9 and 
10. The data at 140 GeV demonstrate a sooth behaviour and Fig.9 is quite conclusive: the 
present model (x 2 = 0.5) provides a much better agreement with the data than the EL- 
model (x 2 = 2.6) or the strong factorization (x 2 = 4.5). Note that the EL-model with n=2 
could not fit the observed mass dependence of the suppression even with a larger k. The 
point is that in the EL-model, the suppression decreases with produced mass because of the 
factor {Qo/Q) n in (4) and this contradicts the trend of the data in Fig.9. The data at 286 
GeV have large error bars in the most important region (x^) 1 / 2 > 0.5 and all theoretical 
curves have x 2 < 1-5 at this energy. However, it is our impression that our results are in a 
reasonable accord with the data at 286 GeV. 

Let us summarize the results of this section. The Figs. 6,7,8 and 10 make no definite 
preference between the models. But the Figs. 5 and 9 allow us to conclude that the present 
model describes the DY production by pions much better than the EL-model. These two 
figures also present an evidence that the strong factorization is violated in nuclei. 

V. DY PRODUCTION BY PROTONS 

As in the case of the DY production by pions, the DY production by protons is dominated 
by the quark-antiquark annihilation subprocess. The difference between these two reactions 
is that protons probe mainly the antiquark content of a target, because protons have no 
valence antiquarks. Therefore this reaction is sensitive to a possible contribution of excess 
mesons in nuclei. The cross section for producing dilepton pairs in proton-nucleus collisions 
is given by Eq.(12). We have used the free nucleon structure functions as q{{x) and q{{ x ) 
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for the following reason. The nuclear modification of the antiquark distribution is unknown, 



but from the electroproduction data[15| it follows that the EMC-like nuclear effects should 



be unimportant in the region X2 < 0.3, measured in the current experiment with 800 GeV 



protons |fL4] . The nuclear shadowing at small x^ is not included in structure functions for 
the same reason as in the previous sections. Thus, in this reaction we will study the net 
effect of the initial state absorption represented by the factor F DY . In the previous workpif. 
we have considered the same reaction with the constant parton absorption cross section. In 
that paper, we have concluded that the excess pion contribution with the average excess 
pion number per nucleon n n m 0.07 is consistent with the data. We will now show that the 
X\— dependence of the absorption cross sections, introduced in the present paper, drastically 
changes the role of initial state absorption in this reaction and affects the conclusion of 



Rcf.[31 



We used the free nucleon structure functions from Ref.||27||. We have integrated over 
the same mass region as in the experiment ||14||. The x 2 — dependence of the ratio R for the 
iron and tungsten targets is shown in Fig. 11 by the solid and dashed lines, respectively. An 
important result is that the suppression is concentrated mainly at 22 < 0.05. The origin 
and the location of this depletion are exactly the same as in the case of the charmonium 
production by protons at 800 GeV (the solid line in Fig. 2). This observation is supported 
by the fact that the calculated R from Fig. 11 is in accord with the data (x 2 = 0.7). The 
A-dependence of the suppression at low X2 is also well reproduced by our model. In this 
situation, there is no room for any noticeable contribution of excess pions, at least at Xi > 
0.04. If our present model is correct, than the upper limit for the momentum fraction carried 
by excess pions is about 1%. This conclusion is qualitatively consistent with the results of 
Ref . [0 that the structure functions of off- mass-shell pions are significantly depleted. In this 
case the excess gluons, considered in Sec. Ill, are another possible candidate to carry the 



missing parton momentum in nuclei |p0 



The xp— dependence of R for the iron is shown in Fig. 12. Though the qualitative trend 
of the data is more or less reproduced by the theoretical curve, the suppression at xp > 0.5 
is overestimated by our model. In that region, the values X2 < 0.04 have the noticeable 
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contribution. We may assume that in the region x 2 < 0.04, not covered by the data points 
in Fig. 11, either the nuclear structure functions are enhanced due to some excess particles 
or our model overestimates the suppression. In any case, our model leaves no room for the 
parton recombination effects in this reaction, as well as in the charmonium production. 
VI. DISCUSSION 

We have assumed a new type of initial state effects, the x±— dependent absorption of 
fast initial partons. As we argued in Sec. I and in Ref . |i~0H , such effects in nuclei are not a 
priori excluded by the factorization theorem forecasts. Using only one physically motivated 
empirical function for the absorption cross sections, we have qualitatively explained nuclear 
effects in several reactions with fast protons and pions. We claim that in the currently 
discussed experiments on charmonium production, the x— dependence of nuclear suppression 
can be accounted for by the initial state interaction we assumed. In this case, the final state 
effects in the region xf > can be fitted by the constant cross section o~f. With the excess 
gluon contribution included, the present model can explain also the energy and projectile 
dependence of the charmonium suppression. The dilepton production suppression is also 
qualitatively explained, though our model sometimes overestimates the suppression. This 
possible discrepancy can be removed by adjusting the parameters, for example by reducing 
a max I a max- The present results are obtained using a minimum of arbitrary assumptions 
and allow us to conclude that this model is able to self-consistently describe nuclear effects 
in different hadropro duct ion reactions. The T— and J/ip— production at low and negative 
xp— demands the incorporation of comover interaction and nuclear modification of gluon 
distribution and will be considered somewhere else. 

If the present model is correct, then the following qualitative conclusions can be made. 
1) The strong factorization is violated by the initial state interactions in hadroproduction 
on nuclei. The weak factorization takes place instead. 2) The parton recombination has 
no significant contribution to charmonium suppression. 3) There is an indication of the 
presence of excess gluons in nuclei at 0.06 < x 2 < 0.15, carrying about 3% of the total 
momentum. 4) The energy-loss model cannot self-consistently describe the quarkonium and 
dilepton production if the mass dependence of the energy loss is given by l/Q 2 . 5) There is no 
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indication of the excess pion contribution in dilepton production in proton-nucleus collisions. 
The fraction of the nuclear momentum, lost by partons due to nuclear binding, can be carried 



by the excess gluons. This point illustrates the difference between the binding correction |3C] 



and pionic models [33] for the EMC-effect: the former is not based on the presence of excess 
pions in nuclei. 

The above conclusions do not mean that the soft initial state interactions, assumed by 
the EL-model, take no place in nuclei. Instead, these interactions may be responsible for the 
observed pr — dependence of hadroproduction on nuclei. But the role of these interactions 
in the x— dependence of nuclear suppression can be minor, in agreement with the limitation 
(3). The validity of the present model can be tested, e.g., by remeasuring with a better 
accuracy the DY production by pions, where we predict several specific pecularities. We 
believe that the results of this paper can be useful for an analysis of nuclear effects in hadro- 
and electroproduction, even if our present model will be ruled out by future studies. 
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Figure captions. 

FIG.l. The function sm 6 (x7r/2), representing the x— dependence of parton-nucleon ab- 
sorption cross sections. 

FIG. 2. The nuclear effectiviness a versus 22 for charmonium production in proton-nucleus 
collisions at the energies 800 GeV (solid line) and 200 GeV (short-dashed line). The latter 
with excess gluons included is represented by the long-dashed line. The parton-recombination 
model prediction is shown by the dotted line. The data at 800 GeV (diamonds) and 200 
GeV (crosses) are from Refs. @ and U, respectively. 

FIG. 3. The same as in Fig.2 versus xp (without the parton-recombination model predic- 
tion). 

FIG. 4. a versus xp for the charmonium production in negative pion-nucleus collisions 
at 200 GeV (long-dashed line). The same with excess gluons is shown by the solid line, a 
for 200 GeV protons without excess gluons is shown by the short-dashed line. The data 
(diamonds for pions and crosses for protons) at 200 GeV are from Ref. fl[] 

FIG. 5. Ratio of the dilepton yield for tungsten to proton versus x<i in pion-nucleus 
collisions at 140 GeV: the solid line is the present model result, the long-dashed line is the 
result of the energy- loss model |J and the short-dashed line is the strong factorization result. 
The data are from Ref. fi~3H . 

FIG.6. The same as in Fig.5 at 286 GeV. 

FIG. 7. The same as in Fig.5 versus x\. 

FIG. 8. The same as in Fig.5 versus X\ at 286 GeV. 

FIG. 9. The same as in Fig.5 as a function of dimensionless dilepton mass. 

FIG. 10. The same as in Fig.5 as a function of dimensionless dilepton mass at 286 GeV. 

FIG. 11. Ratio of the dilepton yields for protons at 800 GeV for iron (solid line) and 
tungsten (dashed line) as a function of x%- The data for iron (diamonds) and tungsten 
(crosses) are from Ref. |14[ . 

FIG. 12. Ratio of the dilepton yields for protons at 800 GeV for iron versus xp. The data 
are from Ref. [IHI . 
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